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ABSTRACT 

We investigate whether or not a supermassive black hole (SMBH) with mass > lO^M© can be made in 
the heart of ultraluminous infrared galaxies (ULIGs) during the course of mergers between/among gas 
rich galaxies, (a) If one progenitor galaxy had a seed SMBH with mass of ^ IO^Mq, this seed SMBH 
can grow up to ^ IO^Mq due to efficient Bondi-type gas accretion during the course of merger given a 
gas density in the circumnuclear region of nn ~ 10^ cm~^. (b) Even if there was no progenitor galaxy 
with a seed SMBH, star clusters with compact remnants (neutron stars and/or black holes) produced in 
the circumnuclear starbursts can merge into the merger center within a dynamical time scale of ~ 10^ 
years to form a SMBH with ;> IO^A/q. Note, however, that the contribution of compact remnants 
supplied from hidden star clusters is necessary to lead to the formation of a SMBH. In conclusion, the 
ultraluminous infrared galaxies observed in the local universe can make a SMBH in their center during 
the course of merger either by gas accretion onto a seed SMBH or by dynamical relaxation of compact 
remnants made in the violent circumnuclear starbursts. Therefore, it is quite likely that the ULIGs will 
finally evolve to optically luminous quasars as suggested by Sanders et al. 

Subject headings: quasars: general - galaxies: starburst - galaxies: star clusters - galaxies: interactions 
- infrared: galaxies 



1. INTRODUCTION 

Since the discovery of ultraluminous infrared galaxies 
(Soifer et al. 1984; Wright et al. 1984; see for a review 
Sanders & Mirabel 1996), these galaxies (hereafter ULIGs) 
have often been considered as possible precursors of opti- 
cally bright quasars (Sanders et al. 1988a, 1988b; Norman 
& Scoville 1988). This argument is based on the following 
observational properties of ULIGs. (a) Their bolometric 
luminosities amount to ~ IO^^Lq, being comparable to 
those of quasars (Sanders et al. 1988a). (b) Their lu- 
minosity function is similar to that of quasars in the lo- 
cal universe (Soifer et al. 1987; Sanders et al. 1988b). 
(c) All the ULIGs are galaxy mergers or heavily interact- 
ing galaxies (Sanders et al. 1988a; Lawrence et al. 1989; 
Leech et al. 1994; Clements et al. 1996). Morpholog- 
ical evidence for galaxy mergers has also been obtained 
for a number of optically selected quasars in the local 
universe although the majority are giant elliptical or gi- 
ant elliptical-like galaxies (McLeod & Rieke 1994a, 1994b; 
Disney et al. 1995; Bahcall et al. 1997; McLure et al. 
1998). However, infrared-selected quasars tend to reside 
in morphologically disturbed hosts (e.g., Hutchings & Neff 
1988; Boyce et al. 1996; Baker & Clements 1997). If 
most giant elliptical galaxies were formed by major merg- 
ers between/among disk galaxies (Toomre 1977; Barnes 
1989; Ebisuzaki, Makino, & Okumura 1991), it is possi- 
ble that the majority of quasar hosts are major merger 
remnants, (d) ULIGs in later merger phases tend to have 
active galactic nuclei (AGN) on the average (Sanders et al. 
1988b; Majewski et al. 1993; Borne et al. 1997). Although 



all the above properties suggest an evolutionary link from 
ULIGs to optically bright quasars, its plausibility is still 
in question. 

It is generally considered that the quasars are powered 
by the central engine of active galactic nuclei (AGN); i.e., 
disk-gas accretion onto a supermassive black hole (SMBH) 
and masses of SMBHs in quasar nuclei are estimated to be 
M, ;> lO^M© (e.g., Rees 1984; Blandford 1990). There- 
fore, if an evolutionary link exists between ULIGs and 
quasars, we have to explain either the presence or the for- 
mation of SMBHs with mass higher than ^ 10* Af© in the 
heart of ULIGs. This issue was already discussed by Nor- 
man & Scoville (1988). They investigated the fate of a 
coeval, massive-star cluster of 4 x lO^M© within the cen- 
tral 10 pc region (see also Weedman 1983) and found that 
a SMBH can be formed in the heart of ULIGs. However, 
recent high-resolution optical and near-infrared images of 
a number of ULIGs using the Hubble Space Telescope have 
shown that the intense star forming regions are scattered 
in circumnuclear regions up to ~ a few kpc from the nu- 
cleus (Shaya et al. 1994; Scoville et al. 1998; Surace et al. 
1998). Therefore, it still seems uncertain whether or not 
a SMBH with ^ lO^M© can be made during the course of 
merger evolution in ULIGs. In this Letter, we investigate 
this issue taking actual observational properties of ULIGs 
into account. 

2. FORMATION OF QUASAR NUCLEI IN THE HEART OF 

ULIGS 

Morphological features of ULIGs suggest that most 
ULIGs come from mergers between or among galaxies 
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(Sanders et al. 1988a; Taniguchi, & Shioya 1998). An- 
other important property of the ULIGs is that they are 
very gas rich; e.g., ~ lO^^M© (Sanders et al. 1988a; 
Scoville et al. 1991; Downes & Solomon 1998). Therefore, 
the progenitor galaxies of ULIGs are rich in gas such as 
giant spiral galaxies. Since intense starbursts are observed 
in many ULIGs, the most probable formation mechanism 
of SMBHs is the collapse of compact remnants of mas- 
sive stars (Weedman 1983; Norman & Scoville 1988). An- 
other important issue is whether or not progenitor galaxies 
had SMBHs originally in their nuclei. Although masses of 
SMBHs in nearly spiral galaxies (i.e., progenitor galax- 
ies of ULIGs) are of the order of 10^~^Mq at most (e.g., 
Kormendy et al. 1998 and references therein), these seed 
SMBHs could grow due to gas accretion in circumnuclear 
dense-gas regions during the course of the merger. There- 
fore, we consider two cases; 1) at least one progenitor had 
a SMBH with M. ~ lO^-'^M©, and 2) no progenitor had 
a SMBH. 

2.1. Mergers between/ among Nucleated Galaxies 

If a progenitor galaxy had a SMBH in its nucleus, this 
seed SMBH could grow in mass during the course of merger 
because the central region of a ULIG is very gas-rich. Such 
gas accretion may be efficient within the central 1 kpc re- 
gion because the gas density is quite high in the central 
region; i.e., ~ lO^""* cm~^ (Scoville, Yun, & Bryant 1997; 
Downes & Solomon 1998). We here consider the classical 
Bondi-type (Bondi 1952) gas accretion onto the SMBH. 
This gas accretion rate is given by 

M = 2 TT mn nu rl v^, (1) 

where mn, ^^h, ?'a, and Vc are the mass of a hydrogen atom, 
the number density of the hydrogen atom, the accretion 
radius defined as ra = GM,v~^ (M, is the mass of the 
seed SMBH), and the effective relative velocity between 
the seed SMBH and the ambient gas, respectively. A typ- 
ical dynamical mass of the central 1 kpc region is of the 
order of ^ 1O^M0. Suppose that the SMBH with mass 
of IQ^Mq is sinking toward the dynamical center of the 
merger. Its orbital velocity is Worb — 67M^^^gr]~^^^ km 

s~^, where Mnuc,9 is the dynamical mass of the central 
1 kpc region of the merger and ri is the radius in units 
of 1 kpc. The crossing time of the SMBH is Tcross ~ lO'^ 
years. Therefore, the merging time scale is estimated to be 
Tmorgor ^ lOTcross ^ 10® ycars (c.g., Barncs 1989). Adopt- 
ing an average gas density in the circumnuclear regions of 
ULIGs to be nn ~ 10^ cm~^, we obtain the accreting gas 
mass during the course of merger; 

Mace = M Tmerger ~ 3 X 10^M.%nH,3t^e"lOO^merger,8 Mq 

(2) 

where M.^e is the mass of SMBH in units of lO^M©, nH,3 
is the average gas density in units of 10"^ cm~^, Wo.ioo is the 
orbital velocity with respect to the ambient gas in units of 
100 km s~^, and Tmerger,8 is the merger time scale in units 
of 10® years. This estimate implies that the seed SMBH 
cannot grow up to M, ^ 10®Mq if the seed SMBH is less 
massive than l(f Mq. Since the ULIGs come from merg- 
ers between two galaxies or among several galaxies, their 
progenitor galaxies should be very giant spiral galaxies in 



order to pile up molecular gas up to ^ IQ Mq in their 
central regions (e.g., Sanders et al. 1988a). In fact, nearby 
spiral galaxies such as M31 and NGC 4258 have SMBHs 
with a few IQ'^Mq [e.g., Kormendy et al. (1998) and ref- 
erences therein; see also Miyoshi et al. (1995) for the case 
of NGC 4258]. Therefore, it seems quite likely that the 
seed SMBH may be more massive than that adopted in 
the above estimate. If the seed SMBH is more massive 
than a few lO'^M©, it could grow up to ^ IQ^Mq. Al- 
though we have no knowledge about the seed SMBHs in 
the progenitors, our estimates given here suggest that the 
gas accretion in the dense gas clouds onto the seed SMBH 
can lead to the formation of a quasar nucleus in the heart 
of ULIGs. 



2.2. Mergers between/ among Non-nucleated Galaxies 

Next we consider a case where there is no seed SMBH in 
the progenitor galaxies of ULIGs. In this case, a possible 
way to form quasar nuclei in ULIGs is to pile up the cir- 
ciimnuclear star clusters of compact remnants of massive 
stars; black holes and neutron stars, each of which has a 
few Mq at most. This issue was pioneeringly discussed by 
Weedman (1983). Using both starburst models by Gehrz, 
Sramek, & Weedman (1983) and optical spectroscopic ob- 
servations (i.e.. Ha luminosity), he suggested that star- 
burst galaxies with L(Ha) ^ 10^^ erg s~^ could produce 
compact starburst remnants up to ~ 10'^ Mq. However, 
since the Ha luminosity of starburst galaxies is dominated 
by the most massive stars in the starburst region, it seems 
hard to estimate the total mass of compact remnants solely 
using L(Ha). Thus, the estimate by Weedman (1983) may 
provide a rough upper limit for the compact remnant mass 
in the starburst galaxies. Norman & Scoville (1988) also 
discussed the formation of quasar nuclei in ULIGs. How- 
ever, their assumption (a coeval, massive-star cluster of 
4 X 1O^M0 within the central 10 pc region) turns out 
to be unlikely because the recent high-resolution optical 
and near-infrared imaging by the Hubble Space Telescope 
of the ULIGs have shown that blue star clusters are dis- 
tributed in the circumnuclear regions up to r ~ a few kpc 
(Shayaetal. 1994; Surace et al. 1998; Scoville et al. 1998). 
Although the central star cluster associated with the west- 
ern nucleus of Arp 220 is very luminous and its mass is 
estimated to be > 10^ Mq, typical masses of the circumnu- 
clear star clusters are of the order of Md ~ 1O''M0 at most 
(Shaya et al. 1994; Scoville et al. 1998; Shioya, Taniguchi, 
& Trentham 1998; see also Taniguchi, Trentham, & Sh- 
ioya 1998). Although some star clusters may be hidden 
because of heavy extinction (Scoville et al. 1991; Genzel 
et al. 1998), we firstly investigate whether or not the star 
clusters found both in the optical and in the near infrared 
(NIR) (Shaya et al. 1994; Scoville et al. 1998) will be re- 
sponsible for the formation of a SMBH with M, ^ 10® M©. 

Shaya et al. (1994) discussed the fate of circumnuclear 
star clusters; since these clusters will lose their kinetic en- 
ergy to individual stars during random encounters (i.e., 
dynamical friction) , they will sink toward the merger cen- 
ter within ~ 10^ years. Here, from a viewpoint of dynam- 
ical relaxation of the star clusters, we examine whether or 
not compact remnants formed in the circumnuclear star- 
forming clusters can make a SMBH with M, ^ IO^Mq. 
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For simplicity, we consider a case where ten eircuninu- 
clear star-forming clusters, each of which has a total stel- 
lar mass of 10^ Mq, are distributed within r ~ a few kpc. 
It is known that stars with m* > 8Mq produce compact 
remnants. We estimate how many such massive stars are 
formed in each cluster. 

We assume that stars are formed with a Salpeter-like 
initial mass function (IMF); 



(3) 



where TO is the stellar mass in units of Mq and /? is a 
normalization constant determined by the relation 



pm, 
J mi 



(j){m)dm = 1, 



which leads to 



1 - (m;/m„)'*-i ■ 



(4) 



(5) 



The number of stars with a mass range mi < m* < m2 is 
estimated as 



N{mi < m 



pm 

* < m2) = / 



TO 



dm. (6) 



Using P in equation (5), we re- write equation (6) as 



N{mi < TO, < 1712) 



(TO]^^ — TO2^) stars Mq^. 



(7) 

There are three free parameters; the power index (fi), and 
the upper and lower mass limits of the IMF (to„ and mi). 
Since stars with to* > 8Mq produce compact remnants, 
we set TOi = 8Mq and TO2 — to„ > mi. In Table 1, 
we give results for some possible combinations of the pa- 
rameters. Although fi = 1.35 is the canonical value for 
stars in the solar neighborhood, there are some lines of 
evidence that massive stars are more preferentially formed 
in such violent star- forming regions; i.e., a top-heavy ini- 
tial mass function with to; ^ IMq (e.g., Goldader et al. 
1997 and references therein). Therefore, we adopt a case 
of = 0.35, to„ = 3OM0 and mi = IMq. In this case, 
there are ~ 4 x 10^ massive stars with to* > 8M0 in each 
cluster. Each compact remnant has a mass from a few Mq 
(for neutron stars) to several Mq (for black holes). There- 
fore, the total mass of compact remnants in each cluster is 
~ 10^ Mq. Such ten clusters will be relaxed dynamically 
with a time scale of 



Ta,n - nirC-y^rTM-,'/' ~ 9.4xl0V^gM-y^ (years), 

where Vci^i is the typical size of a circumnuclear region with 
ten star clusters in units of 1 kpc and M^^r is the mass 
of the clusters in units of IO^Mq. Therefore, we expect 
that a SMBH with mass of ~ lO'^ Mq will be made ~ 10^ 
years after the onset of circumnuclear starbursts in ULIGs. 
This mass is smaller than that necessary for quasar nuclei 
(i.e., ^ IO^Mq). Each cluster of compact remnants would 
be able to gain its mass by gas accretion as discussed in 
section 2.1. However, the accreted mass is estimated to 



be ^ 3 X IO^Mq for ten clusters in total, being still less 
massive than 10^ Mq. 

Here it should be again remembered that all star clusters 
in the central region of ULIGs cannot be observed in both 
the optical and the NIR because inferred extinction toward 
the nuclei of ULIGs is very large; e.g., ^ 50 mag (Gcnzel et 
al. 1998; Scovillc ct al. 1991). Therefore, it is quite likely 
that the majority of nuclear star clusters in ULIGs are hid- 
den by a large amount of gas and dust. Recently, Shioya 
et al. (1998) analyzed the optical-NIR spectral energy dis- 
tributions of nuclear star clusters in Arp 220. They found 
that these clusters are more massive systematically (i.e., 
^ IO^Mq) than circumnuclear ones but can account only 
for about one-seventh of the total bolometric luminosity of 
Arp 220. Although OH megamaser sources found in the 
central region of Arp 220 may provide possible evidence 
for hidden AGN (Diamond et al. 1989; Lonsdale et al. 
1998), the recent mid-infrared spectroscopy of a sample of 
ULIGs has shown that the majority of the ULIGs such as 
Arp 220 are powered by nuclear starbursts (Genzel et al. 
1998; Lutz et al. 1998). Therefore, it is strongly suggested 
that some hidden star clusters should be responsible for 
the remaining (^ 6/7) bolometric luminosity of Arp 220. 
Since compact remnants produced in the hidden clusters 
will join to form a SMBH, it is expected that a SMBH 
with M, ^ IO^Mq will be formed in the heart of ULIGs. 



3. DISCUSSION 

We have shown that a SMBH with mass ^ lO^M© 
can be made in a ULIG during the course of merger be- 
tween/among gas rich galaxies regardless of the presence 
of a seed SMBH in progenitor galaxies; i.e., (a) if one pro- 
genitor galaxy had a seed SMBH with mass of ~ 10 ''Mq, 
this seed SMBH can grow up to ^ IO^Mq because of effi- 
cient Bondi-type gas accretion during the course of merger 
given a gas density in the circumnuclear region of rifj ^ 10"^ 
cm~^, and (b) even if there was no progenitor galaxy with 
a seed SMBH, the star clusters of compact remnants made 
in the circumnuclear starbursts can merge into the merger 
center within a dynamical time scale of ^ 10^ years to form 
a SMBH with ^ IO^Mq. Note, however, that the con- 
tribution of compact remnants supplied from hidden star 
clusters is necessary to explain the formation of SMBHs. 
In conclusion, the ultraluminous infrared galaxies observed 
in the local universe can make a SMBH in their center dur- 
ing the course of merger either by the gas accretion onto a 
seed SMBH or by the dynamical relaxation of star clusters 
of compact remnants made in the violent circumnuclear 
starbursts. 

The presence of a SMBH with M, ^ 10** Mq is a cru- 
cially important necessary condition for the occurrence 
of quasar activity. In the local universe, the masses 
of SMBHs in giant spiral galaxies (e.g., our Milky Way 
Galaxy, M31, NGC 1068, and NGC 4258) are as low as 
M, ^ 10*^^^Mq (Kormendy et al. 1998 and references 
therein). Therefore, it is suggested that isolated, typi- 
cal spiral galaxies cannot harbor quasar nuclei. However, 
mergers between/among gas-rich galaxies can cause effi- 
cient gas fueling toward the nuclear regions of the merging 
systems and then trigger intense starbursts either as a re- 
sult of the piling of a lot of gas (Mihos & Hernquist 1994) 
or by the dynamical effect of SMBH binaries (Taniguchi 
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& Wada 1996; Taniguchi & Shioya 1998). Furthermore, as 
demonstrated in the present work, these mergers provide 
a possible way to form SMBHs with M, ^ IO^Mq. In this 
respect, it is quite likely that the ULIGs will finally evolve 
to optically luminous quasars as suggested by Sanders et 
al. (1988a, 1988b). 

Finally it is worthwhile noting that some elliptical galax- 
ies could be formed by galaxy mergers (e.g., Toomre 1977; 
Barnes 1989; Ebisuzaki et al. 1991). We also note that 
some elliptical galaxies (e.g., M87, NGC 3115, NGC 3377, 
NGC 4261, and so on) have SMBHs with M, ^ IO^Mq 
(Kormendy et al. 1998 and references therein). Actually, 
investigating the physical conditions of ULIGs (i.e., mass 
density and velocity dispersion), Kormendy and Sanders 
(1992) found evidence that ULIGs are elliptical galaxies 
forming by merger-induced dissipative collapse (see also 
Wright et al. 1990; Baker & Clements 1997). Therefore, 
we suggest that almost all SMBHs with M, ^ lO^M© in 
the local universe were made by galaxy mergers. 
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Table 1 

The number of massive stars per unit mass (Mq) 





mi 


ruu 




N{m^ > 8Mq) 
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0.0278 


0.35 
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30 


0.080 


0.0409 


1.35 
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0.437 


0.0189 


1.35 
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30 


0.503 


0.0187 



